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I, INTRODUCTION 


Mie oetay 2S a2 Statistical investigation using 
Empirical Orthogonal Function (EOF) analysis and cross 
correlation analysis to examine the interrelationships 
between anomalous wind forcing and anomalous sea surface 
temperatures (SST) over the North Pacific Ocean from 20°N 
homeo N and 110°W to 130°E for the period January 1969 
monmomel December 1978. Particular attention is given to 
investigating the effect of surface wind components having 
time scales of ten days or less. The primary hypothesis 
to be tested in this research is that the sea surface 
meiig@enature responds to atmospheric wind forcing, as mani- 
fested in storms, and that the amount of the response is 
Peooeretonal to the amount of wind forcing. No other surface 
temperature changing processes are included in this work. 
The concept for examining storminess comes from the strong 
relationship that exists between synoptic storms and air-sea 
exchange [Simpson (1969), Fissel, Pond and Miyake (1977), 
Puoberry and Camp (1978)]. in addition, Davis (1976), who 
examined the statistical relationship between monthly sea- 
level pressure anomalies and SST anomalies, suggested that 
anomalous atmospheric storminess might have a strong corre- 


Phone co anomalous Sol. 
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The approach and expectations of this study were that: 
(1) EOF analysis would be used to extract the essential 
features from the complex space-time patterns of wind and 
popeecatd, 1,€5, it would extract the "Signal" from "noi/se"’ 
mimeeme data, (2) correlation analysis would be used to 
examine relationships between the variables once the noise 
was removed, and (3) anomalous storminess values, computed 
Mimene LOrm Of friction velocity cubed and wind stress curl, 
would be significantly correlated with anomalous SST cooling. 

The atmospheric data for this study were a special 
series of six-hourly surface wind analyses prepared by 
Fleet Numerical Oceanographic Center (FNOC). Friction 
Meme@city cubed (co) ands the Verereal eomponene Of the curl 
of the wind stress (curl, T) were calculated because of the 
eeaolashed relationship between u,> and the mechanical 
mixing of the sea by the overlying atmosphere and the well- 
eeeaolished relationship between curl ,t and Ekman convergence 
and divergence. Monthly mean anomalies of ee and curl t 
Pememcompited USing: (a) high-pass filtered (periods less 
than ten days) wind components only, (b) a combination of 
high- and low-pass (periods greater than 10 days) filtered 
mundmecomponents, and {c) the unfiltered (total) wind com- 
POoMemts, The partitioning of wind parameters into these 
three parts was designed to identify the effect of storms, 


Poulimineisolation and in the presence of the existing mean 


de 








(low-pass) flow. Thus, six fields were created to represent 
the anomalous wind influence upon the ocean surface. EOF 
analysis of all data fields, including the anomalous SST 
foamwes, Was employed to extract the significant and almost 
Meese -tree “geophysical patterns." The resulting anomalous 
wind parameters were then correlated with the anomalous SST 
meeaeat the time lags of minus three through plus three 
months to investigate possible temporal relationships. 

The first topic considered below is data set description 
and analysis. It is followed by a discussion of EOF analy- 
Seoeand adn explanation of the truncation scheme used to reduce 
noise. Then some observations are given on the lagged corre- 
acrons of ey and SST, and curl t ideo Le esta tile tea 
Significance of the correlations is carefully examined and 
iieespretations are made in light of that significance. 
Pinally, some conclusions are drawn about the results 


produced. 
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Il. DATA SET DESCRIPTION AND ANALYSIS 


A. ACQUISITION AND ANOMALY COMPUTATION 
1. Sea Surface Temperature 

Monthly sea surface temperature (SST) data for the 
period January 1969 through December 1978 were obtained from 
Siewsctipps Institution of Oceanography Climate Research 
Group through the generosity of Jerome Namias and Robert 
Born. The geographical orientation of these data as delin- 
eated in Figure 1 was a nine by twenty-five array utilizing 
the 5° square grid used by Namias' group which extends from 
Mieco oy) N and 110°W to 130°E. Due to missing data during 
meme -month record at some of the 225 grid points, only 
166 grid points were retained for analysis. The small 
squares and circles in Figure 1 show these 166 points. At 
these points, only 79 discrete values were missing out of the 
potas I95920 possible. These missing values were filled by 
temporally interpolating between the preceding and following 
monthly values. All SST values were converted from Fahrenheit 
@omccmebearade and then processed to compute anomalies. The 
Some@ratton Of monthly anomalies from ali data sets in this 
aucoemveceaccomplished by finding the ditference between the 
observed monthly value and the 10-year average value found 


at the same grid point and month of year as the observation. 
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Peeeiret1on Velocity and Surface Stress 

iia chetOotmeee: Sol, Six anomalous surface wind 
Parameters were investigated. Three involved re MMe Re wu. 
Peete fr1et10n Velocity and three involved curl t where T 
is the surface stress. These were obtained from Haney, 
Peoomeana Heise (1980), which in turn was based in part on 
Beemeariien thesis by Risch (1980). They used @ special 
series of six-hourly surface wind analyses over the North 
Pacific Ocean prepared by Fleet Numerical Oceanographic 
Center (FNOC) covering the period January 1969 through 
Peeemoer 1978 to analyze for synoptic storm activity. 
Temporal variance of the surface wind components at each 
grid point was "high-pass" filtered, thereby extracting all 
temporal variability having time scales less than ten days. 
iar wchires oe and thee curl + data sets were then calcu- 
lated using, respectively, (a) the total (unfiltered) wind 
components; (b) the high-pass wind components only; and 
(c) a combination of high- and low-pass wind components 
which represent the effects of storms (high-pass) in the 
presence of the existing mean (low-pass) field. See Haney, 
Risch and Heise (1980) for details. Having calculated the 
above wind parameters at every six-hourly observation over 
the entire 10-year period, monthly means, monthly climatolo- 
gies over 10 years, and finally, monthly anomalies were 


constructed. The monthly anomalies, covering the same 


15 








10-year period as the SST anomalies, for the six wind- 
related data setsS were then interpolated from the FNOC 
Perthern Hemisphere Polar Stereographic grid to the same 


latitude-longitude grid points as the SST anomaly data. 


B. DETRENDING AND STANDARDIZATION 
The seven data sets representing monthly anomalies of 


. 


u, Fm Senn (botn from total wind components), vu, Fi, 


d 


curl trl (both from high-pass wind components only), eee. 
eur tr? (both from high- and low-pass wind components) 
Mmcmooat were then detrended and standardized. The detrend- 


ing equation was of the form 


where te is the detrended value at the mee Pinemon Sem am lom, 
Wee Wemede oOricinal observation, a iS the slope, b is the 
intercept and Cc. = (Ci een 1 pee ee COme sion elm 

the time interval. This detrending was performed to remove 
mieiinear trends in the anomalies. 


Standardization of the seven data sets was completed 


by using an equation of the form 
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where X 1S the standardized value, x is the detrended value 
See ene anomaly at a particular grid point, X is the time 
eames valle Of x, and 5 15 the standard deviation of x - x. 
The use of standardized variables assures that each variable 
at each grid point is of equal importance in determining 
the form of the field for comparison to other variable 
fields. Additionally, standardized anomaly fields of sea 
Surface temperature and wind parameters correspond more 
Seocely to the climatological classifications of normal, 
above normal, much above normal, etc., than do the anomaly 
Peebas themselves [Kutzbach, 1967]. Standardization is 
medmemned in Order to utilize the EOP method of extracting 
Pomenae = trom noise," sn Poa CccmeNCt war Loenerie SUDSSQUemne 
correlation analysis. 

The standard deviation maps of all seven data fields are 
@nciluded in Figures 2 through 4. Figure 2 shows that the 
Semeechssonr activity for the three ihe alsicel  Sisie cy eelebingts 
Borwacd ot Japan at about 45°N and 165°E. The largest 
@aemiecs OL the curl fields in Figure 3 appear to be located 
aeemees Neand 170°E. Finally, the standard deviation pattern 


eaeenewoot data 1S indicated in Figure 4 with the most active 


mieavale east of Japan, 45°N and 150°E. 


ie, 
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IfI. CALCULATIONS PERFORMED 


A. EMPIRICAL ORTHOGONAL FUNCTION ANALYSIS 

ietcerscscarch used empirical orthogonal function analysis 
to reduce a large number of highly correlated variables to 
a small number of independent variables. These are linear 
combinations of the original ones and contain a large per- 
wemeape OL the total variability. This technique has been 
used extensively in meteorology and oceanography [see, e.g., 
Peemipden, 1967; Stidd, 1967; Davis, 1976; Weare, 1977; 
pope ho7 2) and, unlike most other orthogonal representations, 
pmlomcechinique does not require a predetermined function, 
mremeas the Fourler components, Sine and coSine, but depends 
upon the spatial and temporal interrelationships within the 
data being analyzed. 

ime zeaen (1967) provides a lucid outline of the mathe- 
matical procedure necessary to define the orthogonal 
milgelons and their coefficients. A covariance (or correla- 
tion) matrix must be formed from the data set. To maximize 
the variance explained by the new functions, one must solve 
the system of linear equations associated with this co- 
Variance matrix for a set of eigenvalues. The eigenvalue 
memmesents the fraction of the total variance in the original 
irae ladsexplained by 1tS Corresponding eigenfunction. A 


function is found by substituting a given eigenvalue into 


Zi 





the system of equations and solving with the aid of an 
additional normalization condition. Just as the sines and 
cosines of Fourier analysis are space-orthogonal, the eigen- 
functions produced by EOF analysis are also Sspace-orthogonal. 
PesOctated with each function is a (time) coefficient which, 
when combined with the function and the result summed over 
eleeunetions, Provides an exact reproduction of the 


Seeoemat £ield. That is, 


2 
mix, t) = S A,,(t)EOF (x), 


where EOF iPmeic ClAcnitunction. ain an 1s the time coefficient 
mea ted to the mth CIgeival ue. ile cemeOct PleTenes. Wwiten 


are time-orthogonal [Kutzbach, 1967], are obtained from the 


: 
functions and the original anomaly field. The specific 
formulations of the mathematics described here are detailed 
In the Appendix. 

Figure 5 shows the cumulative percentage of variance 
explained versus eigenvalue mode number for the seven sets 
of data examined in this study. For all the seven data 
mic@iswmmtne CcUMUlative percent variance tends to increase 
Eocapmenmtrecally with the mode number. Mode number one 15 
mictmomerrtial pattern which, in a least squares sense, 


explains the most variance possible (averaged over all grid 


Pomm@ieseanad time). That variance 1s between 9 and 19% of the 
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Eeerieeveriance depending on which of the seven fields is in 
feeceoumon. it 15 Significant that the cumulative variance 
explained by the SST and ne modes rises more rapidly than 
that of the curl t modes. This means that the major part of 
ficot and u,~ fields Gan be explained by a few Spatial 
modes while the curl variance is distributed more equally 
among a greater number of modes. The next section describes 
the method of determining, for each data variable, how many 


memes tO retain for analysis. 


Pee NCALTTON METHOD AND SIGNIFICANCE LEVEL 
PomereusSsinvestigators [Beale et BO er adda land 
eee loo9; Farmer, 19/71; Dyer, 1975; Rinne and Jadrvenoja, 
1979] have suggested numerous methods for determining the 
Momm@esOr truncation, i.e., finding the sufficient number of 
mempenemes Of the series, given by the empirical orthogonal 
functions, to be retained for analysis. The choice of method 
depends largely on the purpose for which the empirical 
Smeeteconal function eXpansion 1S to be used. In this study, 
eeeetca)! Orthogonal function analysis 1S used to reduce 
Bemsce 's ihe term noise is defined here as that variability 
Seecamaommess confined to small space scales, purely local 
Pee rceend Observational errors. Recent numerical experi- 
fomiesmiyerrersendorfer and Barnett (1977) indicate that the 
classical asymptotic theory of random eigenvalues is inade- 


quate to provide significance tests when the number P of 
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Spatial data points is in the range of 10 to 100 (normally 
encountered in geophysical practice) and when the historical 
meee size N iS in the similar range of 10 to 100. In the 
present study P = 166 and N = 120. Their experiments revealed 
that the asymptotic theory takes hold only for N on the order 
of 1000 or more for the above P values. After numerous 
experiments, Preisendorfer and Barnett have developed an 
emomeeieal procedure for empirical orthogonal function trunca- 
tion in the case of small sample sizes based on analysis of 
Various arrays filled by Gaussian random varlables of unit 
variance, zero mean. The procedure provides, for a given P 
and N, the significant confidence band for the randomly 
generated eigenvalues, thus providing the objective method 
Somraentitvying the "noise’' in a given data set. Tim Barnett 
generously computed and provided, for the specific N and P 
(eiose im tiis study, the 5 and 95% confidence limits of 
eigenvalues tor modes one to twenty-five. Figure 6 demon- 
Strates the use of this truncation process for SST elgen- 
values. The figure shows eigenvalues plotted against mode 
number on a semi-logarithmic grid. The first eigenvalue is 
always the largest, since the first EOF explains the most 
Variance. The two solid curved lines labeled 95% and 5% 
Represent the indicated range of the eigenvalues obtained 
Dolemaicet mrom one hundred repeated calculations using 166 


by 120 arrays of random variables of unit variance. Thus, 
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out of the 100 repeated experiments, 95% of the eigenvalues 
Zemwichnin the rather narrow band indicated by the two 
SUitieope with 54 above the top curve and 5% below the bettom 
curve. The e's represent the eigenvalues of the SST anomaly 
meld, also standardized and, therefore, having unit variance. 
Notice that the arrays of random numbers have the variance 
spread rather evenly throughout the various modes while the 
Poeittcid has much of its variance in the first ten or twenty 
Sime Or the total of 120 modes. According to Preisendorfer 

and Barnett (1977), the method of extracting the "geophysical 
Signal'' from the SST data is to keep those empirical orthogonal 
function modes with eigenvalues lying above the 95% band 
Seeaeed by the random data fields. Im this case, only 5% 
@emene retained field is likely to be attributable to random 
Wemeece sin the case of the SST anomaly data, the first 11 
modes were retained. This truncation procedure was applied 
memeach data field with the results shown in Table I. Table I 
shows, for each data field, the number of empirical orthogonal 
function modes which were retained and the cumulative percent 
of the total variance which those modes explain. 

Newewtnat the data used in this study nave a rather large 
amount of noise. The SST and ee anomaly fields are about 
ioewotenal and 25% "noise," while the curl anomaly fields 
Pence o2 ton65% signal."  Therégis more “noise"™ in the 


curl field because it is a spatial derivative of a quantity 
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(stress) which contains small-scale errors and large-scale 
geophysical variation. 

Figure 7 is a comparison between the SST values at an 
individual grid point before and after EOF truncation. 

Curve (A) represents the detrended and standardized SST 
values from grid point number 84, located one-third of the 
way between Hawaii and the Aleutian Islands. Curve (B) 
Mepresents the resultant values derived from the truncated 
memeecrtes tor the same grid point. It should be noted that 
mepoceits) rerlected in (A) by the curve's rapid small-scale 
meriations. It is readily apparent that curve (B) is 

a damped form of curve (A). A considerable amount of the 
small-scale fluctuations have been eliminated with curve (B) 
Still retaining the large-scale trends of curve (A). 

In addition to the method described above, various other 
ways have been suggested for determining the point of trunca- 
tion of the series given in EOFs. Rinne and Jdrvenoja (1979) 
have offered some of their own methods and summarized many 
others. These methods included: retaining eigenvalues which 
Eeemain more than one percent of the total variance, 
retaining eigenvalues that are much larger than others, 
mmlieeineed statistical criterion based on a Bartlett's test 
to determine the values to be saved, and, finally, keeping all 
the largest eigenvalues which together account for 99% of the 


moetiariance, Another truncation test mentioned in this 
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summarization 1S one which uses a diagram drafted similar to 
Figure 6 in this paper. The procedure truncates the EOF 
series at the largest eigenvalue at which the eigenvalue 
(versus mode number) curve can be approximated with a straight 
mene. in Figure 6 the straight line approximation would 

Start best at mode number twenty-four and thus, according to 
this method, the largest twenty-three eigenvalues would be 
retained. 

All of the above methods would have retained more modes 
than the Preisendorfer and Barnett scheme which was used in 
Mio otudy. Ihe choice of this scheme can best be explained 
by uSing a schematic of Rinne and Jadrvenoja (1979) shown in 
meeure §. 

iitemwiele box represents the total variance and this 
@emotsts Of real environmental information (areas 1 and 2) 
emamnoise (areas 3 and 4), separated by the interior vertical 
line. EOF series truncation controls the placement of the 
Sloping broken line which divides the box into an explained 
(areas 1 and 3) and unexplained (areas 2 and 4) portion. The 
mumoered areas are now identifiable as: (1) real information 
pomictimedsby EOFS, (2) real information EOFs fail to explain, 
Mopmemioice explained by EOFs, and (4) noise that EOFs fail to 
eeelladme (A) represents the situation when a relatively 
small number of EOFs are retained while (B) occurs when a 


relatively large number is saved. The Preisendorfer and 


ou 








acu eo, 


B 


BOMeImterpretation Of the total variance where 
A is the case of saving a small number of terms 
and B is the case of retaining a large number 
of terms in the EOF series expansion. 
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Barnett (1977) method used in this study is closely represented 
by (A). Very little noise is retained at the "expense" of 
explaining a somewhat reduced amount of total variance. In 

the following section, possible relationships between the 


seven truncated anomaly data sets are examined. 


oe 





ee, EN ee tes 


Pe LAG CORRELATIONS BETWEEN oe veo: 

Figures 9 to 15 are maps of lagged correlations between 
oe and SST anomalies (detrended, standardized, and truncated) 
over the North Pacific Ocean. The map series steps through 
lags of minus three to plus three months, in one-month 
increments, with Figure 9 showing the minus-three-month case 
and Figure 15 ending the series with a plus-three-month lag 


- 


between tee and SST. In the map series, SST is lagged 


relative to a. Pits ttcoeMepev ieee lag Of Minwsescuree 


months shows the correlation between the fe seen) ocr Clic 
Seer mtiat occurred three months before the mee Tine ac 
opposite situation, the map lagged plus three months shows the 
correlation between the we field and the SST field that 
mecimeod tiree months later. Each figure contains three 
maps, which are, from top to bottom, correlations between 
SST and (1) fal computed from high-pass wind components 
Omiy, (2) 2 computed from a combination of high- and 
low-pass wind components, and (3) u,° Fm computed from 
unfiltered {total) wind components, respectively. Finally, 
each map is shaded gray in areas of negative correlation 
values, 


iieamereipation that anomalous stoOrminess values are 


foe wedetomanomalous Surface cooling, It was expected that, 
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Figure 9, Correlation maps of the a anomaly data 
fields with SST at a lag of minus three 
months. The last closed isoline is labeled 
in some regions having the largest correlations. 
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Same as Fig. 9 except lagged plus one month. 
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Figure 14. Same as Fig. 9 except lagged plus two months. 
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Figure 15. Same as Fig. 9 except lagged plus three months. 
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ror lags greater than or equal to zero, maps would be mostly 
covered by negative correlation values. Examination of 

1 1S u.> Pel acomoOeeiiawomtevedls thdtea sitenif icant change 
occurs in the zero and plus-one-month lag maps. There 

1S a visible increase in the gray shading, i.e., negative 
correlations, at these lags. Table II delineates some map 
Pememetcrs that cam be used to quantify this increase in 
the negative correlation area when SST lags ee Almost 
without fail there are significant changes in going from 
M@eeiinis-ene-month jag to the plus one. For instance, 
the sum of the negative correlations (RNEG) for Pe 82 

goes from -10.81 at minus-one-month lag to -30.19 at plus- 
One-month lag. The number of grid points with negative 
correlations (NNEG) jumps from 99 at minus one to 145 at 
zero and then up to 147 at plus-one-month lag for Boe 
The significance of this jump can be appreciated by con- 
Eeaerinoead fandom correlation situation. In such a case, 
foe ieepoints would be expected to have positive correla- 
tions and 83 would be expected to have negative correlations, 
Peeeuecmsctamaadrd deViation Of 6.5 grid points, This is 
Paaterieewerue Only if the grid points are independent of 
each other. This means that 68% of the time the number of 
negative points would lie between 76.5 and 89.5 and 95% of 


the time the number of negative points would lie between 70 
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mmamgO., thus, the increase of negative grid points from 99 
momla7 1S Of notable significance. 

BetOucwcomeinuing beyond this point, the issue of 
feoniticance must be examined further. An issue of signifi- 
cance was raised and resolved earlier in this study when 
@Bensiacring where to truncate the EOF series to ensure 
maximum signal and minimum noise. Thus, the correlations 
presented above are those associated with the geophysical 


7 


fitmietuations in SST and Wee Cer ello me cre enak ) rand, are sao 
feeurbutable to error ‘'noise."' The paragraph preceding this 
Mmemmoted the significance of the changes that occurred in 
the correlation between negative and positive lags and pro- 
vided mathematical evidence for this significance. The 

Meer oniticance to be examined is related to the existence 
of artificial correlation due to the interplay between the 
femeth Of record examined and the dominant time scales of 
waeiatiton Of the processes being studied. 

Pomaetyrst Approximation to the magnitude of this 
Peeolote artificial correlation, Six grid points were chosen 
femeeles im Fig. 1) and their temporal history of SST, com- 
Muiseadeatter EOF truncation, Was examined for repetitive 
@vents. (See Figures 16 through 21.) SST data were chosen 


Somectimate the magnitude of the possible spurious correla- 


tions because it has relatively longer time scales of varia 


Mibieercempared to the atmospheric fields. Table III shows 
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Biemvalwes Obtained from these calculations. The number of 
independent events, n, was estimated by averaging the number 
of peaks with the number of troughs in the time series. 
Dividing 120 months by n provided the average event length. 
iitesrtormulia for Rcrit, the critical correlation value, was 


derived from Davis (1976) in the following form, 


po 


Rerit = ey 





3] 


where n is again the number of independent events during 
Meemaata record length. At any grid point, a correlation is 
Seemordered significantly different from zero only if it 
Peeeeeoas herit. the importance of these values of Rcrit 
Peeeeeoe discussed in a detailed examination of all the corre- 
faeton maps in a tollowing paragraph. 

Sciremammano (1979) has described a method of presenting 
@eerelations which conveniently and accurately incorporate 
Mamormoninheance Of the correlation relative to the level 
Peretticial correlation. Thé technique consists of 
Momtairzime (dividing) the raw correlations by their large- 
Migestandard error (O.). iio i@mcams non tne Bante tt 


pOoumula (Bartlett, 1978, p. 352] as 


CAG) Cyy CiAt)/N (4-2) 


a2 





where . and oo are the discrete autocorrelation functions 
Gf the variables x and y, respectively. It can be seen that 
(4-2) represents the ratio of a measure of the integral 
time scale of the variables x and y to the record length N. 
Thus, the normalization will compensate for the interplay 
of the dominant time scales of the input processes and the 
record length and will yield significance levels. Considera- 
tion of this interplay can be crucial when comparing meteor- 
ological and oceanic variables which, in general, have 
feecerent time scales of variability. In the present study, 
Such a presentation of the correlations would essentially 
Sive a more quantitative estimate of Rcrit than given above, 
@mad it would allow Rcrit to vary over the North Pacific, 
reflecting the fact that the dominant time scales of varia- 
eemelety Of SST and ee ayaa Det ago em Orme mr UMens oleae) (eneste 
Calculations were not made for this thesis and would offer 
AO pper tunity in a subsequent study to increase the value 
mthe correlation maps. For the purpose of this thesis, 
pmectrelation, r, is considered Significantly ditferent 
meomezero if |r} > 0.20. 

iemibteht of the previous comments, further review or the 


tee lag correlation maps and Table II will now reveal several 
Paictonabeinteresting details. All maps retlect a large 
geographic scale of correlation compared to the total area 


of interest. The timing of the maximum magnitude of negative 


a 





area (NNEG) and maximum total negative correlation (RNEG) 
Varies among the three data sets. Table II shows that 
Merl has its maximum area of negative correlations at the 
zero month lag, while u,F2 and u, Fm Wave Cie ly Massa weds 
established at a plus-one-month lag. The maximum map sum 
Sema negative correlations occurs at a plus-One-month lag 
OT heel and ava? and at a zero month lag for u, Fm, The 
values of RNEG and NNEG attained by the eee Gidea at aa 
of plus one month exceeded all other data maps. Thus, it 
appears the u,> Vomlrbe set Ne re wOCe radeon taat anoma |ous 
eeermess values, as reflected by ieee omen relanea Ge 
ememalous Surface cooling with a time scale of two to six 
weeks, 

B@piceermeerestimes tO note that if only correlations of 
Pepeeure value greater than 0.20 are considered significant, 
then most of the trends indicated above remain in effect. 
For diem and Eee tie pbk Omeaence DOcweeI ve  COrre lation 
@Gueads are not significant, while approximately half of the 
meet ve correlations (at lag of +] month in Figure 13) are, 
The result appears to be, after considering all the possible 
elemtetcance checks provided previously and the map details 
delineated above, that at a plus-one-month lag there is a 
iiueeemid-oceanic area where ae aaa and ae eee ieee das 
anomalous storminess variables produce a cooling of the sea 


SUrface, The mechanism responsible for this cooling most 


24 





imeeely ts ethe vertical mixing (entrainment) of cold water 


from below the mixed-layer. 


B. LARGE LAG CORRELATIONS BETWEEN ceed AND SST 


Somat experiment to lend additional credibility to the 
large shift in correlations observed at the plus-one-month 
lag, an analysis was conducted by correlating the SST field 


5 


[acted wind parameter, specifically u,° Fl, that occurred 


, 
fey, 21, 24, 29, and 36 months later. It was expected 
Miates the correlations obtained from this analysis would be 
non-significant and offer by their scale and magnitude some 
Support for the previously described results. Figure 22 
momeuins thescorrelation maps, from top to bottom, of minus 
memminis 24, and Minus 36 months’ lag and Figure 23 contains 
@ieweerrelation maps, from top to bottom, of minus 17, 
Minus 21, and minus 29 months' lag. Examination of these 
Gyo ,rigures reveals that, in comparison to eer lagged plus 
one month, all maps have much smaller areas of negative 
eomretations, although some areas are significant (greater 
mutes 20) by the earlier test. 

apie IV which delineates the values of RNEG, NNEG, 
RPOS, and RSUM for u,°Fl at the trial lags examined indi- 
cates that the RNEG and NNEG values remain much below those 
values observed for cea at zero and plus-one-month lag. 


The average RNEG is -7.67 and the average NNEG is 76 for 


lags of minus one to minus 36 months. Both values lie 
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Bigumes 22. Correlation maps ox u, Fl and SST Anomaly 
fields for lags of minus 12, 24 and 36 
months. The last closed isoline is labeled 


in some regions having the largest correlations. 
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Figure 235. Same as Fig. 22 except for lags of minus 
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considerably below the values observed for Teed ates = 


one-month lag, where RNEG was -24.96 and NNEG was 131. 


C. LAG CORRELATIONS BETWEEN curl NID oo 

mteures 24 to 30 are maps of lagged correlations 
between curl t and SST over the same area as the previously 
viewed correlation maps of ore and Sook the layout (rime 
maps 1S exactly the same as before with the maps progressing 
through lags of minus three to plus three months, in one- 
month increments with Figure 24 showing the minus-three-month 
case and Figure 30 ending the series with a plus-three-month 
lag between curl t and SST. In this map series, as in the 
Meemmneus series, Sot 1s lagged relative to curl 7. ase 
the Map with a lag of plus one month shows the correlation 
between the curl, t wield -aneethe woolen cld stmt mOCCULT Ca some 
Meme tater, Each figure contains three maps, which are, 
magmecop tO bottom, correlations between SST and (1) curl TFl 
computed from high-pass wind components only, (2) curl tF2 
computed from a combination of high- and low-pass components, 
and (3) curl tFm coliputed from uUnfmiltered = total) wind com- 
mememes. respectively. As before, seach map 1s shaded gray 
Pmimemeaseort mecative correlation values. 

Simla expectations existed for this set of correlation 
maps like those indicated before for ie Thats. Sainec 
poomalousmevalues of Wind stress curl indicate anomalous 


Ekman pumping, it was expected that for lags greater than or 
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eroememe4. Correlation maps of the curl anomaly data 
fields with SST at a lag of minus three 
months. The last closed isoline is labeled 
aesOne Ge eLonsmhavence the Fanreestucorrelations. 
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Figure 25. Same as Fig. 24 except lagged minus 
two months. 
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Piece 20, Same as Fig. 24 except laoced minus 
one Month, 
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Piemnenc/, oame as Fig. 24 e€xeept nowlac between 
data sets. 
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Figure 28. Same as Fig. 24 except lagged plus one month. 
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Ewes). came as Fig. 24 eXNeGCDimmrgcced plus’ three 
months. 
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equal to zero, the maps would be mostly covered by negative 
aeemcliatton values, Examination of Figures 24 through 30 
reveals that the correlations between auGde Fidos eoeceun On 
a smaller geographic scale than the correlations between u.> 
waeeool. lhere is an inerease in the amount of gray shaded 
area as the maps are viewed from minus three to plus three 
months lag with a maximum occurring at the zero or plus-one- 
month lag. This increasing gray shading trend is not as 
Serons as that obsrved in the ae series and examination of 
Table V shows that the statistics used to quantify the 
correlations do not approach those values obtained from 
tape II, the nae ie eouretemuc we neldomem, sreveals. fidt aire 
max number of negative correlations (NNEG) is attained by 

ini curl «Fm field with a plus-one-month lag. The maximum 
magnitude of totaling all negative correlations (RNEG) is 
Seeerved tO OCcUur at a plus-one-month lag in the curl 1F2 
field. Remembering that the random correlation situation 
would have a standard deviation of 6.5 grid points, both 
curl tF2 and curl tFm undergo shifts of approximately two 
Standard deviations or more when the lag shifts from minus 
emiemeonzere and from zero to plus one month lag. Thus 

there is a significant shift away from the possible random 
correlation case. When only the statistically significant 


Gomnetacions (iri > .20) are considered, much smaller 


Moree eirercourelation areas result im Figure 28. These 
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Pitciericant areas are Small and Scattered. in the gulag EZ 
meee a negative Correlation area persists off the West 
Coast of Canada with a maximum size and strength at plus one 
Maem ceneral, the Significance of the relationship 
between SST and curl t appears to be somewhat less than that 


~ 


between SST and ah 
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Vee BQUNG IU SON 


The primary hypothesis that was tested in this study 
was that the sea surface temperature (SST) responds to 
atmospheric wind forcing, as manifested in storms, and that 
the amount of response is proportional to the amount of wind 
Forcing. When the atmospheric forcing is characterized as 
omeamomalous storminess and calculated as a and curl 7, 
mamas a negative correlation with anomalous SST over most 
of the North Pacific at a plus-one-month lag. It appears 
that ee P-wwOkew> eros lvwamec lated topos! anomalies than 1s 
curl t. aise reread! Mieing May bewa Stronger influence 
than Ekman divergence and convergence in surface cooling. 
Peewee thtsedaiterence could be caused by difficulties 
iiemeTan im COMputing curl [Saunders, 1976]. 

Pmemner conclusion concerns the e€ompensation required 
m@mvalidate tne magnitudes of the map correlations. The 
memetively snort length of record (120 months} necessitates 
pPeeolcula@tiron Of the artificial correlation and a ssignifi- 
fimo reduetron 1m the computed correlations. Any increase 
ieee cecerd length will reduce the level of artificial 
@eeeolation 210r this type of study [see (4-1)}], but to make 
Geeemomiricant increase in the total validity, a doubling or 


Cimtewotineno: the total record lengthews mecessary. 
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The present study examined the relationships between 
PiecmeemOspnere and Ocean that exists throughout the year. 

It is now well known [Davis, 1978] that many ocean- 
atmospheric relationships are strongly seasonally dependent. 
imaes  pOSSible that if the data o£ this study were stratified 
mieowsecasens, Statistically significant results could be 
Smeained. the problem of artificial correlation would 
mm@enrease because of fewer independent events [n in Table [Il 
@mid {4-1)] available for correlation, but the actual corre- 
Meerons COUld poOSsiDlyeincrease enormously. For example, the 
memeatiOn ©: the Winter anomalous storminess with anomalous 
sea surface temperature cooling has the potential to be much 
more Significantly correlated than the non-seasonal relations 
Peetoped in this research, 

Whenever one endeavors in a study of this sort, more 
questions are often generated than answers found. One 
jmmiimnine question is, What can the wind alone tell about 
the future state of the ocean? That is, perhaps vee and 
Gurl >t are sutficiently uncorrelated that they Ee gerner 
[meee wmere Sieniticantly correlated with SST. It is 
recommended therefore that a further study correlating SST 
With a combined form of tte and curl, t Memacine in Ordeal © 
provide a more complete answer to this question. This will 


give a more complete estimate of the total importance of 


Giemiinasmmebelative to other exieécts, in changing the sol. 


ro 





eee ND 


Methods for representing arrays of environmental 
variables by empirical orthogonal functions have been 
formulated in a number of different ways. [See, for example, 
meomeema., 1956; Kutzbach, 1967; Stidd, 1967; Davis, 1976. ] 

Mest Of the information presented in this appendix is based 
Meom the work of Kutzbach (1967). 
itAm— mom ScOnwononEevoeLor WiLGM represents the 


n 
7 ie wee MP ObsServation Of P variables and let A 


1 


Bema PF by IN Matrix, where the n*" column is the observation 
mee tor ee ieecmcOad le tS tOcetermimesa vector € Which has the 
highest degree of similarity to all the AD observation 
Meeters Simultaneously. Kutzbach (1967) explains that if 

the similarity 1S-measured witn the squared and normalized 
inner product between the observation vectors Re and the 


determined vector e, the problem is to maximize 


go (e'A)(e'A)' e'AAte 
1* 1 = a a — i eae 
aoe ee Wee nee 
after averaging across all AL: (new op gune scenotes the trans- 
pose, This is equivalent to maximizing 
Soa (1-1) 


(a 





Siegecet tO the condition 


? Gl =Z) 


Pete 6B iS an m by m Symmetric matrix whose Bea element 


Mmemwervyen by 


B= N+ (A'A) , es) 


and mis the larger of P and N. 
SeiemMeme.  NUtZbacn (1907) continues by pointing out 
aieewtaesMaximization of (1-1) subject to (1-2) leads to 


the equation 

Be =ed , (1-4) 
The vector e and the parameter » can now be recognized as 
an eigenvector and associated eigenvalue of the symmetric 
normalized matrix B. It can be shown that the e, are 
orthogonal and the ds are real and positive. Therefore, 


(1-4) can be rewritten as 


BE = EL (1-5) 


LS 





where E is the m by m matrix of eigenvectors and L is the 


m by m diagonal matrix of eigenvalues. Note that 


E’E = 1. (1-6) 


Bomeie following, it 1s assumed that the elements of L and the 


columns of E have been arranged such that e, is associated 


1 
with ho: the largest eigenvalue of B, C5, 1s associated with 
re cue mMemcemlareest ecicenVvalite ot Be and so forth, Combining 


Wes) ancoew(1=5)2 and noting from (1-6) that for orthogonal 


iererieces the transpose 1S equal to the inverse, one obtains 


E'AA'E = LN Cle) 


defining 


=e = foe (1-8) 


where C is an m by m matrix, it follows that 


> 
i 


EC Gls) 


and 


N. (i-10) 
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Examination of (1-10) reveals that the observation vector 
AL can be expressed as a linear combination of the m eigen- 
Meetors. It all m eigenvectors are used, then the observa- 
Eom vector ae 1S eNACtie mop meduced, 

Pmnoste applications Gf emparical orthogonal function 
Bemilon-random data, it 15 found that a large portion of the 
Momtance Can DE accounted for by retaining only the first 
few largest eigenvalues and related eigenvectors. Thus, 
retaining only the first few terms of (1-10) approximates 
mere, OobServation vector. If A” is defined as the approxi- 


h 
mation to the n°” observation Ve G GO weet 
Ces. ares oe en: PC Ce 


This allows truncation of noise by removing those terms with 
small variances. These smaller terms are most likely to 
consist almost entirely of noise as shown by the work of 
Peedsendgrier and Barnett (1977). 
prcemmmuemano ll cation Of the preceding 15 delineated 

Peomemisime oo! data, in an outline detailing the array 
berzecmanage'bi-560 IMSL routines. Note: Only the trunca- 
MiOieeritts (the Value of Pp”) Vit icdmrniOieeuine data Sets and 


Pe elgomorta points with N = 120 observations. 


3 





ARRAY 


A 


eiz0 &K 


166) 


A'A (166 x 120)(120 x 166) 


i. 


AC 


wlo6 x 


Ploo. x 


(166 x 


(166 x 


(e720 Xx 


POO 120 


166) 


166) 


iy 


2} 


elz0) X= 1o6) 


NAME 


Anomalous 
Datea oer 


Transpose 
Proqduee 


Normalized 
SY imine tec 
Ma tric 


Eigenvalues 


Eigenfactors 


rune aved 
EROChiieetOm.> 


Time 
COCERUGrenes 


Approximated 
Hata et 
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